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Directional solidification of lamellar eutectic structures 
submitted to uniaxial stress is investigated. In the spirit of an 
approximation first used by Jackson and Hunt, we calculate 
the stress tensor for a two-dimensional crystal with triangu- 
lar surface, using a Fourier expansion of the Airy function. 
The effect of the resulting change in chemical potential is 
introduced into the standard model for directional solidifica- 
tion. This calculation is motivated by an observation, made 
recently [I. Cantat, K. Kassner, C. Misbah, and H. Miiller- 
Krumbhaar, Phys. Rev. E, in press], that the thermal gra- 
dient produces similar effects as a strong gravitational field 
in the case of dilute-alloy solidification. Therefore, the cou- 
pling between the Grinfeld and the Mullins-Sekerka instabil- 
ities becomes strong, as the critical wavelength of the former 
instability gets reduced to a value close to that of the latter. 
Analogously, in the case of eutectics, the characteristic length 
scale of the Grinfeld instability should be reduced to a size 
not extremely far from typical lamellar spacings. Following 
Jackson and Hunt, we assume the selected wavelength to be 
determined by the minimum undercooling criterion and com- 
pute its shift due to the external stress. In addition, we find 
that in general the volume fraction of the two solid phases is 
changed by uniaxial stress. Implications for experiments on 
eutectics are discussed. 
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FIG. 1. Schematic setup of a directional solidification ex- 
periment. A container with the melt in it is pushed through 
a thermal gradiend G with a velocity V. 



A non-hydrostatically strained solid in contact with its 
melt or vapor can partially relieve its elastic energy by 
producing an undulated interface. This is the cause of 
a morphological instability giving rise to the evolution 
of grooves with a definite spacing under uniaxial stress 
and, possibly, island formation, if the stress is biaxial. 
The instability was first predicted by Asaro and Tiller [|) . 
Experimentally, it has been observed and studied by Torii 
and Balibar ||. Since the independent rediscovery of the 
instability by Grinfeld || , it has often been referred to as 
the Grinfeld or Asaro-Tiller-Grinfeld instability (ATG). 
Important contributions leading to a broad interest in 
the instability are due to Nozieres P|J|. 

In directional solidification (FigjTj), it is known that 
the moving front undergoes, depending on the growth 
velocity, another morphological instability, named after 
Mullins and Sekerka || (MS), where the interface devel- 
ops a cellular structure. Cantat et al. investigated 
the coupling between these two instabilities for dilute al- 
loys. 

They discovered that under favourable circumstances 
a weak uniaxial stress of the order of 1 bar leads to a 
dramatic change in the stability range of the Mullins- 
Sekerka instability. A schematic representation of one 
of the most common liquid-solid equilibrium phase di- 
agrams is displayed in Fig [|. Dilute alloy means that 
the concentration of the minor phase is very small. The 
other situation, in which we are interested here, corre- 
sponds to a composition close to the eutectic one. The 
growing solid then often forms a parallel array of the two 
coexisting phases a and /3 that grow side by side. This 
growth mode is called lamellar eutectic growth. 

A seminal theoretical desription of lamellar eutectics 
has been given by Jackson and Hunt ( JH) Q . Their ba- 
sic idea is the replacement of the diffusion field in the 
liquid phase with that of a planar front. Assuming that 
the a and (3 lamellas have equal average undercoolings 
they were able to obtain an analytic approximation for 
the average undercooling of the interface. They then 
invoked the hypothesis, which has since become known 
as minimum undercooling assumption, that the selected 
wavelength of the pattern leads to the minimum possible 
value of the undercooling (which means that for given 
undercooling the fastest- growing structure is selected) . 



1 




C e C |a C s(3 



FIG. 2. Generic phase diagram of a binary eutectic. T is 
the temperature, c the concentration of one component. The 
regions L, a, and j3 correspond to one-phase equilibrium states 
of the liquid, the solid a, and the solid /3 phases, respectively. 
L+a and L+/3 are regions of two-phase equilibrium between 
the liquid and the solid phases; the actual concentrations of 
the two phases are given by the liquidus and solidus lines 
(full lines) delimiting these regions. c c , c a , and cp denote 
the equilibrium concentrations of the liquid and the two solid 
phases at the triple or eutectic point. The concentrations for 
the undercooled case are also displayed. 



II. MODEL EQUATIONS 

In describing the problem by a macroscopic contin- 
uum model we must introduce fields. These are the tem- 
perature, the concentrations and the stress fields. We 
make some standard simplifying assumptions about the 
properties of the sytem, believed not to affect its essen- 
tial physical features. These simplifications were justified 
elsewhere Juj. For the sake of completeness, we recapit- 
ulate them briefly. The thermal gradient G is assumed 
constant in the frame of reference moving along with the 
growing interface. This means that thermal diffusion is 
much faster than chemical diffusion, that thermal con- 
ductivities of all phases are equal, and that latent heat 
production can be neglected. Thanks to this approxima- 
tion, the motion of the temperature field is completely 
decoupled from that of the concentration field. Temper- 
ature is given by position, which effectively reduces the 
number of fields to be considered by one. We further sup- 
pose the attachment kinetics at the solid-liquid interface 
to be fast on the time scales of all other transport pro- 
cesses. This assumption is legitimate for microscopially 
rough interfaces. We take surface tension to be isotropic. 
In the vicinity of the operationg point in the phase dia- 
gram, the slopes of the liquidus and the solidus line are 



assumed constant. This leads to temperature indepen- 
dent partition coefficients for both phases a and /?. The 
partition coefficients k a m are the ratios of the slopes of 
the liquidus and solidus lines, respectively. In addition, 
we restrict ourselves to the so called one-sided model, i.e., 
we have no diffusion in the solid phases. 

Introducing a dimcnsionless concentration field c — 
(c— c e ) / Ac, where c stands for the physical concentration 
and Ac — cp — c a is the miscibility gap, we can write the 
equation of motion in the laboratory frame (where the 
sample is pushed at constant velocity V along the — z 
direction) 
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In this equation, / = 2D/V is the diffusion length, where 
D is the diffusion constant. One boundary condition for 
the diffusion equation takes into account that the con- 
centration far away from the surface has a constant value 
c oo = (cqc ~ c e )/Ac. In the lateral direction, we assume 
periodic boundary conditions: c{x,z) = c(x + \, z). Mass 
conservation requires boundary conditions for the normal 
derivatives of the concentration fields at the liquid-solid 
interface. This continuity equation reads 
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where 8 = (c e — c a )/Ac is the reduced miscibility gap 
of the a phase and 1 — 8 that of the (3 phase. v n — 
{2D jl + ((x))n z is the normal velocity of the interface 
where the normal points from the solid into the liquid. 
For the stress field we impose mechanical equilibrium, 
ddij/dxj — 0, which means that on the time scale of 
the concentration field, the stress is always relaxed. We 
assume linear elasticity and an isotropic solid, so that 
Hooke's law reads: 
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where <xy are the components of the stress tensor and 
Uij = ^(dui/dxj+duj/dxi) those of the strain tensor (ui 
is the displacement vector). E a {Ep) is Young's modulus 
for the a {(3) phase, v a (vp) the Poisson number. 

The boundary conditions at the solid-liquid interface 
are 



<j nn = nan = -pi 
o nt = nat = , 



(4) 



where n (t) is the normal (tangential) vector at the inter- 
face, and pi is the pressure in the liquid. These conditions 
state that we have no shear at the solid-liquid boundary 
and that the normal component of the stress tensor is 
continuous. That is, we neglect the capillary overpres- 
sure present when the interface is curved. Usually, this 
is a good approximation. 
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There are two further points that have to be taken into 
account. Both result from the requirement of local ther- 
modynamic equilibrium at the interface, due to fast inter- 
face kinetics. The hrst of these is often referred to as 'me- 
chanical' equilibrium condition for the surface tensions of 
the three interfaces meeting at a triple point (although it 
is indeed a condition of thermodynamic equilibrium un- 
der particle exchange, i.e., one of chemical equilibrium). 
The contact angles fla/p (see FigJ|) should obey 



7 Q / sinz9 Q + -ipi sin §p = j a p , 

Jal COS -d a - 7,3; COS <&p = , 



(5) 



where 7^ is the surface tension between the phases i and 
j (and I designates the liquid phase). 
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FIG. 3. Illustration of a lamellar eutectic. The interface 



position is z — C( x )- The pinning angles ■ 



a/0 



are also shown. 



The second condition couples the stress to the concen- 
tration field. It is a modified Gibbs-Thomson equation: 



Ca/P c|i ntcr facc = C/fa"^ + do°^« 
jj a /f3 (Pit - On, 

(e Q = -1 , ep = 1) . 



(6) 



In this equation, £ is the z coordinate of the liquid-solid 
interface and k its curvature, taken positive where the 
solid is convex. 1%/ 18 are the thermal lengths, given by 

= m a /pAc/G, where m a {rap) is the absolute value 
of the slope of the liquidus line describing coexistence of 
phase a (/?) and the liquid. d l = 7izX , c / Lim.iAc are the 
capillary lengths (i = a,/3), where Li is the latent heat 
per unit volume and T c the eutectic temperature. The 
modification is the inclusion of the stress term with 
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Herein, ao is the uniaxial prestress that can be controlled 
in experiments. A detailed derivation of eq. (0) is given 



III. JACKSON-HUNT THEORY FOR A FLAT 
INTERFACE 
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FIG. 4. The flat-interface structure used in the simplest 
Jackson-Hunt approach. 



The first level of approximation in Jackson and Hunt's 
approach consisted in replacing the actual diffusion field 
in (||) with that of a planar lamellar structure sitting at 
the average position of the solidification front. Without 
the stress term, (H) would then become a pair of second- 
order differential equations with boundary conditios fol- 
lowing from (j|) . The solution of these equations with the 
supplementary condition that the two solutions match at 
the triple point gives the interface shape and the volume 
fraction r\ of the a phase. Since these equations are non- 
linear, they cannot easily be solved analytically. Hence 
Jackson and Hunt invoked the condition of equal aver- 
age undercooling of the two solid-liquid interfaces, which 
fixes the free parameter r\ and allows to obtain an ana- 
lytic relation between the average undercooling and the 
wavelength. The second step - solution for the interface 
shape - can then be done numerically, if desired. 

The main modification in our case is that we have an 
additional term in (^|) involving the stress distribution 
at the interface. In the spirit of Jackson and Hunt, we 
compute this expression for a flat interface first. Then the 
problem becomes very similar to JH's original approach 
with the diffusion field replaced by c\i — e a /pH a ^(au — 

(T nn fl(JQ 2 . 

Averaging the diffusion field obtained by solving the 
von-Neumann problem (Q), (Q) for a flat interface, we 
have 
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and the segregation coefficient k has been taken equal in 
the two phases. The averages of the curvature of the a 
and (3 lamellas can be obtained without approximation, 
as they just involve the integration of a derivative, 



(k)q = — r sin?? Q , 
77A 
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To average the stress terms, we must, in principle, solve 
the elastic problem for a flat lamellar structure. Innocent 
as this problem may look, it is not all that trivial. Nev- 
ertheless, the final averaging procedure will turn out to 
be independent of the subtleties that we will now discuss 
briefly. 

At each lamella boundary between the a and (3 phases, 
see Fig. [|, we have, on the one hand, continuity of the 
normal and shear components of the stress tensor (due 
to mechanical equilibrium): 



o xx (x = ) = (7 XX (X = + ) , 

a xz (x = 0~) = a xz {x = + ) , 



(13) 



and the same conditions at x = r/X. On the other hand, 
coherence of the interfaces between lamellas imposes ad- 
ditional conditions, viz. continuity of the displacements 
(up to a constant): 



(14) 



u x (x = ) = u x (x = + ) , 
u z {x = 0-)=u z (x = 0+) , 



with again identical conditions at x — rjX. Equations (\lc 
and (|l4|) and their counterparts at x — rjX constitute two 
boundary conditions at each vertical boundary for the 
stress field in the lamella extending between x = and 
x = r/X. (There are four equations but each of them 
pertains to two lamellas.) The four boundary conditions 
at the two x — const, boundaries of a lamella suffice to 
solve the elastic problem uniquely. Therefore, there is no 
room left for more boundary conditions. But in fact, we 
have, at the boundary towards the liquid 



<7zz(z = C) = ~Vl 

cr xz (z = C) = , 



(15) 



two additional boundary conditions, rendering the prob- 
lem overdetermined. Note that this line of reasoning pre- 
supposes different clastic constants in the solid phases. If 
all elastic coefficients are equal, then the validity of Jl4| ) 
implies that of (|l3|) simply by virtue of Hooke's law (as- 
suming, as usual, that continuous physical functions are 
also continuously differentiable). With different sets of 
elastic constants in the two phases, we have a situation 
similar to that in microstructures discussed by Muller 
JyJ . A solution to the elastic problem need not exist. 
That is, the elastic problem may not have a solution with 



the boundaries fixed to the chosen positions. However, a 
solution to the mathematical problem given all the dis- 
cussed boundary conditions does exist, if we allow the 
lamella boundaries to adjust their shape, i.e., if we con- 
vert the question to a free-boundary problem. The pur- 
pose of the following discussion is then only to establish 
that analytically tractable homogeneous-stress solutions 
exist in particular cases. 

In fact, we do not need general solvability to consider 
a sensible physical problem. Looking for constant-stress 
solutions of (jl4|) together with ( |l5| ) we obtain, setting 
Cxz(C) = Oj the conditions 
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where the superscript indicates the absence of spatial 
variation inside the lamellas and the subscripts a and j3 
distinguish the elastic constants in the two solid phases. 
There are no such subscripts on the stresses and on u° zz 
which are equal in the two phases (in contrast to u xx , 
which may differ). It is evident that for different elastic 
constants in the two materials, ( |l6| ) has a unique solution 
for a xx and u° zzl providing the coefficient determinant 
does not vanish. That is, we just have to choose the right 
value of the prestress o~ xx to ensure the existence of a 
homogeneous solution on which we can base our analysis 
|Q. As long as pi ^ 0, we have a xx ^ — pi, i.e. the 
Grinfcld instability is potentially activated. For pi = 0, 
on the other hand, we can even have a continous set of 
solutions, if we choose the elastic constants such that the 
coefficient determinant vanishes (which is possible even 
for E a ^ Ep, say). 

Given the fact that there is a solution to the elastic 
problem, the calculation of its influence on the Gibbs- 
Thomson equation (|^) becomes very simple. As a xx 
is homogeneous throughout the sample and because of 
o~tt — o~ xx for a planar interface, we simply have {pu — 
o~nn) 2 /o~o = 1- Hence, the averaged stress terms are sim- 
ply H a and H 13 , respectively. 

Inserting this in the Gibbs-Thomson equation, we get 
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where ( ) JH is the average without the stress term. 
Assuming equal average undercoolings in front of both 
phases, we set (() a = ((}p, (because AT = — GQ. As 
has been discussed earlier, this assumption is not nec- 
essary to obtain closed equations [0, but it simplifies 
calculations. We can then write an implicit equation for 
77: 
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(18) 



The last term in this equation is small for small under- 
cooling (implying small Peclet number A/7) and small 
contact angles, so that in this limit an explicit formula 
for i] is available. Using (|l^) in (|l7|), we obtain for the 
averaged undercooling: 



<AT(A)> = (AT) r 



A 



where 



JH 



(V) 



(AT) mln = (AT) 

i 



JH 



G 



l T a l T 



h a + 1 



(19) 

(20) 
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Because on setting d(AT)/dX = the elastic terms dis- 
appear from the equation for A m i n , there seems at first 
glance to be no effect of elasticity on the selected wave- 
length. But that is not true, because rj has changed. 
Expanding A,„i n about T] JH , setting 77 = r\ ,1H + Ary, we 
obtain 



Amin An 



1 P'(T) 



JH\ 



■ (22) 



<" J ">( 1 + a "[-2P(,-, 
do sin $p — do a sin $ a . 
■q JH d o sin dp + (1 - rj JH )d a sintfj 

where A77 w k(l T H^ - l T a H a )/{l T a + l T )- 

The first thing to note is that if the elastic constants 
and the latent heat per volume are equal in the two 
phases, elastic effects do not influence the wavelength 
at minimum undercooling, within the flat-interface ap- 
proximation. This is why we insisted on considering the 
more general case in spite of the complications concern- 
ing the existence of a solution to the elastic problem. 
The logarithmic derivative P'(rf)/P(rj) of the JH func- 
tion vanishes for r\ — ^ and diverges for r\ —> or r\ — > 1, 
allowing for a potentially large effect. However, it stays 
smaller than 50 for 0.04 < 77 < 0.96, which means that 
it does not provide more than an order of magnitude in 
most situations. The second term in the brackets of 
usually is on the order of one. The sign of the effect de- 
pends on the sign of A77, i.e, the relative magnitude of 
the elastic constants in the two phases. 

If we assume that the difference in Young's moduli in 
the two phases is on the order of 10% of their average 
(i.e., (1 - vl)/2Ep - (1 - vl)/2E a « 0.05/£ av ), we find, 
for typical values of the material parameters (T e ~ 400K, 
a freezing range mjAc ~ 10K, L t ~ 10J/cm 3 , k ~ 1, 
E = 10 5 N/cm 2 ) and for 77 - 0.1 that Ar)P'(r))/P(r)) s» 
This gives a relative wavelength 
1 bar and one of 10% for <tq = 100 



2 x 10~ 7 [cm 4 /N 2 ] crj 
change of 10 -5 for <jq 



bar. We therefore conclude that this effect is small in 
ordinary experiments but might be accessible in high- 
pressure setups, where pressures of 100 bar or more could 
be applied. 

The next task is then to see what is the order of magni- 
tude of the influence of deviations of the interface shape 
from planarity. 



IV. JACKSON-HUNT THEORY FOR A 
TRIANGULAR INTERFACE 




FIG. 5. Simplified surface structure 

The simplest non-planar surface structure accessible to 
an analytic approach is a triangular surface (see Fig||). 
To proceed, we will from now on assume that the elastic 
constants are the same in the two phases. 

In the absence of volume forces, the two-dimensional 
stress tensor can be expressed via an Airy stress function 
X- Setting 



d 2 X 
dz 2 



= — 



d 2 X 
dzdx 



we automatically satisfy the condition of mechanical 
equilibrium J"^ daij/dxj = 0. Hooke's law together with 
the assumption of isotropic elastic properties then implies 
that x must obey the biharmonic equation A 2 ^ = 0. We 
split the Airy function according to x(x, z) — x^ z) + 
X^ (x, z), where 
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(24) 
(25) 



K n = 27m/ A, and both terms are solutions to the bihar- 
monic equation separately. 

Equation (24J) corresponds to a homogeneous stress 
state and (p5) describes the deviation therefrom. Once 
we have calculated the coefficients A n , B n we are able 
to compute the stress term in (^|). Inserting our bound- 
ary conditions for the stress field into a representation 
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of a nn and <r„ t in the xz coordinate sytem, we arrive at 
an infinite linear system of equations that in principle 
could be solved for the coefficients. An analytic result 
can be obtained, if the equations are expanded in terms 
of A a m — 2ha/p/^a/0> where X a — rjX and Xp = (l—rf)X 
are the widths of the lamellas and h a (hp) is the height of 
the triangle in the a (/?) phase (Fig. . If the expansion 
is performed up to linear order, one arrives at 



A n = -a a Xe-" m 'A aP n ( V ) , B n = , 

where 

A a p n ( v ) = S n>0 f-A a | + ±A aV 2 + ±Ap(l - v ) 2 



(26) 



[A Q + A^(-1)"-(A C 



A /3 )cos(7r??n)] . (27) 



Note that in ( pq ) we need this definition only for n > 
0, where it simplifies to the second term. Using these 
coefficients in x an d calculating the average of (a t t — 
) 2 , we obtain 



(((T tt — fJ nn ) 2 ) a = (To 2 
((fftt - 0-nn) 2 )/3 = <7 2 

where 
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To be consistent, we have to compute the average of the 
diffusion field for the double trianglular surface as well. 
It turns out that the result can be cast into a form that 
is very similar to the case of a planar interface. All that 
has to be done is to replace the Jackson-Hunt function 
P( V ) by 

I ■s--\ smirrjn 



P(r],A a ,Ap) = P( V ) + -J2 



tt* *■ — ' n 

n—1 



from the formulas. In particular, we assume (Q a ~(()p — 
\{h a — hp). Next, we write down the total average 
undercooling. In minimizing it, we suppose a weak 
A dependence of 77, which yields dP(r],A a ,Ap)/dX = 
-Pi (77, A a ,Af))/\ with A (77, A Q , A p ) = P(rj, A ai Ap) - 
P(rj). We then find that, surprisingly, the result for the 
wavelength does not contain the modified Jackson- Hunt 
function anymore but just the original one: 



E 



sin nrjm 
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and here all integer values, including zero, can appear in 
the superscript of A"^ m (?7). Whereas P(~q) is essentially 
independent of A, the wavelength dependence of 77 being 
weak, P(77,A a , A^) does depend on the wavelength via 
the A dependence of the A a /p. This must be taken into 
account in the minimization procedure when the mini- 
mum undercooling is determined. 

Thus replacing P(n) with P(r), A a ,Ap) in (§) and (g), 
we can proceed in a pretty straightforward manner. First 
we use an assumption analogous to the equal undercool- 
ing assumption to eliminate the term \(coc + 6 + 77 — 1) 
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where 
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1 — 77 \ 77 1 — 77 ' 



(33) 



For comparison with the stress-free case we rewrite this 
as 

A mi „ 2 = A- 2 (77) (l + (^-(l-v)H^(r,) " 
m '" W ^ 2[do Q (l-^)smtf Q +do / Vntf / 3], 

(34) 

where we have taken the Jackson-Hunt result for the 
wavelength at the pertinent value of r\. Of course, there 
is an additional effect (as in Sec. Ill) due to the change 
in the volume fraction under external stress. The latter 
is given by 



A77 = 



It 13 l T a \ 2A 



Ulr p H p -l T a H a ) + ~(hp-h a ) 




P 1 (r ] ,A a ,Ap) 



(35) 



7 p 1 a 



(A— (77)-A^™(77)), (31) 



1 — 77 77 



In order to get an estimate of the order of magnitude 
of clastic effects, we note that for itq « 1 bar and the 
material parameters considered in section III, we have 
fl"^ ps 2 x 10~ 5 . ^(77) is on the order of ten, hence 
0(77) « A, if we take the heights h a ^ of the lamellas to 
be of order A/10. Assuming d^ 13 w 10~ 3 A, we find that 
the second term in (|34|) is on the order of one percent for 
(To = 1 bar, i.e., an appreciable effect may be expected 
for pressures or tensions in excess of 10 bar. 

With the same assumptions, we note that the change 
of 77 induced by elastic effects is on the order of 10~ 4 for 
do = 1 bar and 10 -2 for (To = 10 bar, hence negligible in 
most cases in comparison with the direct effect given by 
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. Of course, this also depends on the size of d\^ n / drj, 
which we have estimated to be small for r\ values not too 
close to or 1, in Sec. III. 

We now consider a few special cases that are especially 
transparent. 

If the lamella structure is symmetric under an ex- 
change of the a and (3 phases, i.e., r) = h and h a = hp, 
then we see immediately from ( |33| ) that Cl(rf) — 0. Terms 
with even n vanish because of the factor sin(7r?7n), terms 
with odd n produce a factor of zero inside the brackets. 
Therefore, application of external stress will not alter the 
wavelength in this case, except possibly via the change 
in r\ induced by (35), which is a much smaller effect. 
Moreover, if we assume the thermal properties of the two 



phases to be the same, i.e. 



T> 



Lp, we have 



H a = H 13 according to (0) (because we took the elastic 
properties of both phases equal from the outset of this 
section). Therefore, we have An = in this case. The 
direct effect on A as described by (|34|) is then absent even 
if h a ^ hp, although there will be a small shift in n, if 
the two phases have different heights. 

Another simplification arises, if we choose all the prop- 
erties of the a and (3 phases to be equal and set A a = 
Ap = A but allow for T) ^ h. In particular, this means 
that we assume the heights of the lamellas to be pro- 
portional to their widths. We can then evaluate f2(r?) 
analytically, 



n( )= 8A y sim^) 



1 + (-1)" - 2C0S(7TO77) 



dx In | sin(7ra)| I , (36) 



8A . , 
= ?yln2 

7T 



and it is easy to show that (2n — > 0. Therefore, 

we have an increase of the wavelength in this case. 

A discussion of the general case is most easily done by 
numerical evaluation of ( p2| ) for a few characteristic sets 
of parameter values and graphical representation of the 
result. This is carried out in Fig. ^. We compare the 
n dependence of the relative change in wavelength for 
A a = Ap, A a = 2Ap and = 2A a . A a is set to 1/10 
and the pressure is 25 bar. The diffusion length is taken 
to be I — 10 2 A and the capillary length d a = 10~ 3 A. The 
contact angles have been chosen as d a /p = arctanA a /^ 
in keeping with the spirit of the triangular approxima- 
tion. It is seen that when there is an asymmetry between 
the lamellas, a decrease of the wavelength can occur, but 
the magnitude of the effect is pretty small if A a w A^. 
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FIG. 6. The change in wavelength A as a function of the 
volume fraction r\ for 25 bar. The thick symmetric curve 
is for A a = Aa = 0.1. The thick asymmetric curve is for 
A a = 0.1 and Aa = 0.2 and the thin curve is for A a = 0.1 
and Aa = 0.05. i? a /a = arctan A Q /a is assumed. 



V. SUMMARY 

To conclude, motivated by the fact that the interaction 
between the Grinfeld and Mullins-Sekerka instabilities is 
strong in directional solidification of dilute alloys [py§, 
we were led to investigate the influence of uniaxial stress 
in directional solidification of lamellar eutectics. 

From the outset, two differences could be expected. 
First, the basic lamellar structure is not determined by 
the MS instability, so direct visibility of an interaction 
with the ATG instability was not likely. Second, since 
the lamellar spacing is typically an order of magnitude 
smaller than cell spacings in dilute alloys, the influence 
of the ATG instability which at typical thermal gradients 
is "resonant" with the MS instability should be expected 
to be weaker in eutectics. 

On the other hand, it is also known that qualitative 
features that are present in dilute alloys, such as parity 
breaking or the appearance of asymmetric cells, invari- 
ably turn up in eutectics, too, albeit often via a different 
mechanism, which is a rather fascinating phenomenon by 
itself. Parity breaking, for example, can be explained by 
two-mode coupling in cellular growth but requires quite 
a different analytic approach in the case of eutectics jl3) . 
More basic features, such as the underlying symmetries, 
are the same in the two cases. 

A similar situation arises here: The mechanism, by 
which stress modifies the properties of the system is en- 
tirely different from that of the dilute-alloy case. There 
it was the coupling to the MS instability, here it is a 
coupling to the asymmetry between the two solid phases. 
Uniaxial stress has a direct effect on the volume fraction 
of the phases, which in general results in a (small) influ- 
ence on the wavelength of the pattern. In addition, it 
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changes the undercooling of the front in a wavelength- 
dependent manner, provided there is a (geometric) dif- 
ference between the a and (3 phases. Both effects were 
calculated to linear order in the deviation A of the front 
shape from planarity. The first effect is present even for 
a planar interface, if the elastic constants of the two solid 
phases differ, and it has been evaluated for that case as 
well. 

As expected, appreciable wavelength changes require 
stresses that exceed those necessary in dilute alloys by an 
order of magnitude. So we do not expect elastic effects 
to strongly affect directional solidification experiments 
with eutectics by accident (which might however happen 
for dilute-alloy experiments). Nevertheless, stresses of 25 
bar or so are not too high to be imposed in a controlled 
experiment which then would allow to test this theory. 

Another point worth mentioning is that the wavelength 
change can be both positive and negative for eutectics 
(and is positive most of the time) whereas we have only 
seen a wavelength decrease with dilute alloys so far (at 
small pulling velocities, the case considered here). This 
makes the effect somewhat less interesting for material 
processing purposes but underlines the basic difference 
in the mechanisms by which stress modifies microstruc- 
tures in the two cases. Large stresses (> fOO bar), how- 
ever, might be used to engineer the volume fraction of 
the phases - if they can be sustained in an appropriate 
experimental setup. 

This work was supported by a 'PROCOPE' grant in 
the framework of a French-German cooperation. 
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